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Abstract: We find spherically symmetric solutions in the modified Horava-Lifshitz gravity proposed 
recently by Bias, Pujolas and Sibiryakov. The non-linear equations of the two derivative action turn 
out to be similar to those stemming from the four-derivative action explored recently. We analyze 
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where the cosmological constant is non-zero. We derive the large distance expansion of solutions 
and show that the power of the standard Newton's law is modified in the presence of a cosmological 
constant. 
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1. Introduction and results 

Hofava provided recently with an alternative view on the quantization of gravitational theories, Q. 
By abandoning Lorentz invariance he was able to allow for a power counting renormalizable class of 
theories. Several aspects of such theories have been analyzed since, and such progress was reviewed 
in where we refer the reader for references and more details. A generic problem of the theory 
proposed by Hofava is the presence of unstable modes as well as strong coupling directions, as it 
usually happens in non-diffeomorphism invariant theories. It was proposed in 0] that the addition 
of spatial derivatives of the lapse function can ameliorate the behavior of the theory, and indeed it 
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was shown that the linearized theory around flat space has proper behavior for a range of the extra 
couplings. 

The purpose of the present paper is to study non-linear aspects of the theories proposed in [[J by 
studying its spherically symmetric solutions. It turns out that the resulting equations for arbitrary 
couplings are much more complicated than those studied (and solved) recently in J|. To proceed we 
will simplify the theory so that it is still tractable while it contains the new ingredients proposed in 

For this we will choose the action to be 

dttfxy/gNlaiKijK* -\K 2 ) + i(',j(" J • + a^a 1 )] , a* = -j— (1.1) 

and we will look for spherically symmetric solutions with zero shift of the form 

ds 2 = - N ( r fdt 2 + ^- + r 2 dn 2 (1.2) 
f( r ) 

Defining the dimensionless variable b = the linearized stability of the theory constraints < y < 2 
or b > 2. 

The nonlinear equations are solved and we find the following types of solutions: 

• In the stable range, b > 2, we find an asymptotically flat solution with positive "mass", which 
extends down to a regular horizon 1 (defined as the solution r* to /(r*) = 0). The asymptotic 
expansion of this solution is 

_ 1 2GM (2GM) 2 (2GMf | 
r 2br 2 Abr 3 

2 2GM 2(GM) 3 /n/ 4x , . 

N 2 = l- — + y- + O(r-) (1.4) 



2GM = i= 7 V ^ r, (1.5) 

This expansion defines the parameter we call "mass" . There is an infinite series of subleading 
terms controlled by the mass and the parameter b. In particular, the subleading behavior in iV 2 
has the same power as the one induced in RS compactifications. There is however an important 
qualitative difference in the sense that in RS the coefficient is linear in the mass while here it 
is cubic in the mass. Therefore, the effects of the new interaction are more prominent for very 
massive objects. As we vary the parameter 2 < b < oo, the ratio r*/2GM in ( Jl .5 ) varies from 
1 (at b — oo where the interactions are not present) to oo at b = 2. 

There is another class of solutions with non-trivial asymptotic behavior 

as r — > oo. For such solutions, the four-dimensional scalar curvature is singular as r — > oo. 



lr This will be a horizon only for regular relativistic matter, see the discussion in 
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In the borderline case b = 2, we find solutions with f(r) arbitrary and 

rN' e 



n - 1 + 7T ' 6 =1 (L7) 

This is in tune with the degeneracy of the spatial derivative action of the scalar mode in this 
case, 0. 

In the unstable range, < b < 2, we also find an asymptotically flat solution with positive 
"mass", which extends down to a regular horizon. Its large distance expansion is similar to 
L~3| ), (1.4). There are no exotic asymptotics in this case. 



• We also analyze the case where a cosmological constant is present. Although we cannot solve 
the non-linear equations in that case, we can solve them for the large distance expansion. We 
classify the large distance asymptotics in that case. We find the regular solutions have an 
expansion that is given by 

f = 2--^— + - t _ ar 2 + ^l_ + ... o = 2-4l (1.8) 

2£ — ai £(3£ — 2ai)(2£ — ax) r ^+2 a\ 

where / is the standard blackness function and C is a constant of integration. We observe that 

q — 2 

the 1/r tail of the solutions without a cosmological constant is now changed to l/r a + 2 . This 
implies that the power of Newton's law is modified in the presence of a cosmological constant. 

Such solutions can be understood as long range tails of matter distributions (stars). The presence 
of a physical singularity in such solutions needs to be further understood. As argued in [0 the 
curvature invariants are not enough to ascertain singular physics. This has been also shown recently 
[|] by analyzing spherically symmetric solutions in the projectable theory. 

A theory with action ( |1 . 1| ) differs from General Relativity already at the two- derivative level. 
Although this difference is at the level of gravitational (and not matter interactions) we expect 
constraints on the coefficient a\ to be mild. We can estimate that the relative change to the perihelion 
precession for Mercury is about — 110~ 8 . Since b > 2 in the stable range, such a constraint is 
comfortably irrelevant due to the non-relativistic suppression. The most stringed constraints on b 
will be set by binary pulsar data and are expected to be in the 10~ 3 range. 



2. The Hofava-Lifshitz gravity theory and its generalizations 

We review here the gravitational theories of the Hofava-Lifshitz type with general couplings. 

The dynamical variables are N,Ni,gij, with (mass) scaling dimension 2 0, except iVj that has 
scaling dimension 2. This is similar to the ADM decomposition of the metric in standard general 
relativity, where the metric is written as 

ds 2 = -N 2 dt 2 + g^dx 1 + N i dt)(dx j + N j dt) , N { = gij N j . (2.1) 

2 Note that this is the canonical scaling guaranteeing power-counting renormalizability once the dimensionfull con- 
stant in front of dt 2 is set to one. 
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The most general power-counting renormalizable action is 

S — Skin + S3 + S re i + S new (2-2) 

S kin = ocj dtd 3 x v ^N[(K ij K^-XK 2 )} (2.3) 

S rel = J dtd 3 x^N [ 1 £ ijk R u V j R l k + CR l] R %J +r]R 2 +iR+(y] , (2.4) 

S 3 = J dtd 3 x^gN [pCijC ij + faRDR + (3 2 R 3 + (3 3 RRijR ij + faRijR^Rh] (2.5) 

S new = J dtd 3 x^/gN [^(aia 1 ) + a 2 (aia 1 ) 2 + a 3 R tj aiaj + a4RVia l + (2.6) 



The extrinsic curvature is defined as 



Kii = ^(9ii-V i N j -V j N i ), , K = g*K ij , K ij = g ih g jl K kl (2.7) 



with covariant derivatives defined with respect to the spatial metric g^. £^ k = the standard 
generally covariant antisymmetric tensor, e 123 is defined to be 1, and other components are obtained 
by antisymmetry. Indices are raised and lowered with the metric g^. Therefore, £^ k = (±l)/y/g. 
The Cotton tensor in fl2.5|) is given by 



C" ----£ ikl V k (Rh-\RVi) (2.8) 



In ( |2.6|) we have defined 

a* = -fi- (2.9) 



and the ellipsis in (|2.6| ) refers to dimension six terms involving aj as well as curvatures. 
The action in (|2.2|) is invariant under a restricted class of diffeomorphisms 

t' = h(t) , a/* = #(*,£). (2.10) 

The transformation of the metric under infinitesimal diffeomorphisms {h ~ 1 + k, h l ~ x % + e 1 ) is 

Sgtj = die k g jk + dje k g ik + e k d k g i:j + kg^ (2.11) 

5N i = d i 6 j N j + d j e j N i + e j g ij + kN i + kN i , 5N = e j djN + kN + kN. (2.12) 

The equations of motion stemming from the action ( |2.2| ) are formidable, even for simple ansatze 
like the spherically symmetric ansatz that will be considered here 3 . We will therefore simplify the 
theory, keeping however some of its main features. We will concretely study the action 

S = J dtd 3 x^N[a(K ij K ij -\K 2 ) + [3QjC ij + £R + a^a^)] (2.13) 



3 The only ansatz that is still tractable is the homogeneous and isotropic cosmological ansatz whose solutions are 
the same as those discussed in M. 
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which contains the kinetic terms, the leading conformally invariant dimension 6 terms, the Einstein 
term necessary for the IR limit and the leading relevant term of the new kind advocated in ||. It 
should be mentioned that the last term breaks Lorentz invariance in a important way and therefore 
we expect its coefficient to be constrained. As it is in the gravitational sector the constraints are 
expected to be mild. 

We will search for static, spherically symmetric solutions with zero shift of the form 



ds 2 



-N(rydif + 



f(r) 



(2.14) 



For such solutions the kinetic and Cotton terms do not contribute to the equations of motion. 



3. The equations of motion 



We insert the ansatz ( p. 14 ) into the equations of motion. The equation stemming from the variation 
of N becomes 

1 



rN' r 2 N' 2 r 2 N" 

2£ (1 - / - rf) - a x (4/ + rf) — + aj-^- - 2aJ- 



N 



N 2 



N 







The equation stemming from varying gu is 

,2 



r 



r N' r 2 N' 2 

4£(/ - 1) + 8£/^f + 2a x f 



N 



while the one stemming from g 2 2 variation is 



-2£r/'-2£(r/' + 2/) 



rN' 



N 2 



2 N n r 2 N I2 

4^f— + 2a 1 f^ r 







(3.1) 



(3.2) 



(3.3) 



w 



a x -2i 



a + b 



N N 

All other equations are trivially satisfied. 

It is convenient at this point to define the following parameters 

For a positive Planck scale we must have £ > 0. Stability of linear fluctuations around flat space 
implies < y < 2 that translates into b > 2 and consequently a < 0. The general relativity equations 
are recovered in the limit a x — > or b — > oo. 

We now proceed to solve ( |3.1|) for iV" and substitute in ( [3. 3D to obtain the equation 



(3.4) 



r 2 N 



12 



rN' 



(a 1 -2O^/ , + 2e 2 (l-/) = 



(3.5) 



N 2 ~~^ J N 

We now combine this with (|3.2| ) to eliminate N and obtain a differential equation purely for /. This 
is valid as long as a ^ 0. We will return to the a = case later. The equation for / reads 

r 2 f 2 + 2arff - Arf + 2a/ 2 - 2(2 + a)f + 4 = (3.6) 

The strategy is to first solve (|3.6| ) and determine /, and then solve (|3.2| ), determine N' /N alge- 
braically in terms of / and then integrate once to obtain N(r). As there are three equations for two 
unknown functions, the solutions of the procedure above must be finally substituted to the original 
equations in order to verify whether there are further constraints. 
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4. Large distance expansions 



Before proceeding with the general solution of ( |3.6|) we will investigate the asymptotic expansion of 
asymptotically flat solutions when they exist. Assuming that at large r / asymptotes to a constant 



(that we normalize to one), we can show that the equation (3J:) admits the following large distance 
expansion 

f = i _2^_{2GMf_{2GMf +(D{r _ i) 
r 26r 2 46r 3 

Ar2 2GM 2(GMf ml 4x . . 

N " = 1 " — + "W" + ° {r ) (4 ' 2) 

where GM is an integration constant that we interpreted as the Newton constant times the mass. 
Note that the corrections to the Newton's law are controlled both by the mass and the parameter b 
that is controlled by the non-standard term in the action (|2.13|). 



5. Finding Solutions 

To proceed further we define 

so that (B.6I) becomes 



We further define 



to rewrite (|5.2|) as 



f = l + - + g(r) (5.1) 
a 



r 2 g' 2 + 2a(rgg' + g 2 ) + 2bg = (5.2) 
g = r ~ a h , r = e u (5.3) 
h 2 + abh 2 + 2be au h = (5.4) 



This equation is similar but not identical to the one solved in . It turns out that a similar technique 
applies here, namely we define 

W = i (5.5) 
h 

so that 

2b e au W 2 2b W 2 

1 + abW 2 ' 9 ~~l + abW 2 { } 

where we have used ( |5.4|) . 

The N function can be obtained from (13.21) that can be rewritten as 



with solution 



rN'\ 2 tN' bf-1 i \ 



rN' b r n „ 

— = --[! + eaW } , e 2 = l (5i 
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Substituting (5.1), (|5.6|) , (|5.8|) into the three equations Q3.1|) -( j3T3D we obtain that that they are solved 
iff W satisfies 

2eW = - (1 + aeW) (1 + abW 2 ) (5.9) 

We note that the equation for e = 1 above is obtained from the equation with e = — 1 by W — > —W. 
As such a sign change does not change /, we restrict ourselves for now on, without loss of generality, 
to the e = — 1 case. 

To proceed further and integrate ( |5.9|) we have to distinguish three cases 



6. Study of the solutions for a < 



r 




w 



Figure 1: Solid lines: r/r$ (vertical axis) is plotted against W (horizontal axis), for a = —1, 6 = 3. There are 
four branches that reach asymptotic infinity atW = - A and W = i — . r reaches zero when W = J= . 

I a l \a\b y\ a \b 

Dashed line: rjr§ (vertical axis) is plotted against W . Thick line: The four- dimensional scalar curvature (in 
units ofr^ 2 ). 

This is the regime where the linearized theory is well behaved. In this case the solution is 



1 + ^\a~\bW 



i+- 



\b+\a\ 



1 + Ial W 



1 - y/\a\b W 



| a | b+ 1 a | 



r 
r 



(6.1) 
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w 



Figure 2: N/Nq (vertical axis) is plotted against W (horizontal axis), for a = —l,b = 3. 



while 



N 1 - ^\a~\bW 
iVo = 



1 / b 

2 V FT 



(6.2) 

1 + y/\a\bW 

In figue |l] the solution W(r) is plotted in inverted form: r/r as a function of W. N/N is plotted 
as a function of W in figure |[ 

Rewriting the solution in terms of g(r) we must distinguish two cases 

• W = 



\a\g-2 



in which case we obtain 



1 + 



\a\g-2 



| a | £>+ 1 a | 



1 _i_ * / M / \ a \d 

1 ~T~ " b V |a| 9 -2 



|o| fl -2 



|a| b+ 1 a| 



r 
r 



(6.3) 



iV 



1 - 



\a\g-2 



1 + 



\a\g-2 



2 V FT 



(6.4) 



The two branches are described in figure |^. Both have a region at infinity and as we will now see 
they have non-standard large distance asymptotics. 

For r — > oo we must have \ag\ — > oo. Expanding around this limit we obtain: 




Figure 3: Left figure: r/r$ (vertical axis) is plotted against g (horizontal axis), for a = — 1,6 = 3 when 
W > 0. The left branch has g < while the right branch has g > A . Asymptotic infinity r — > oo is at 

g ±oo. in the left branch r$ < r < oo u>/iiZe in the right branch \J~^ r < r < oo. The right branch 
corresponds to m < W < oo, while the left branch to < W < ™. Right figure:N /Nq (vertical axis) is 
plotted against g (horizontal axis), for a = — 1, 6 = 3. 




Figure 4: Left figure: r/r$ (vertical axis) is plotted against g (horizontal axis), for a = —1,6 = 3 when 
W < 0. The left branch has g < 0, the center branch has A > g > rr while the right branch has g > A. 

Asymptotic infinity r — > oo is at g = ^ . In Z7ie Ze/t branch < r < ro, in i/ie center branch \J~^ r Q < r < oo 

while in the right branch < r < oo. 27ie Ze/t branch corresponds to A= < < 0, the right branch to 

~~ jaf < ^ < — "ypjk anc ^ ^ e cen ^ er branch to —j^ > W > — oo. Right figure: N/Nq (vertical axis) is plotted 
against g (horizontal axis), for a = — 1,6 = 3, W < 0. 

The metric is therefore not asymptotically flat. Moreover as seen from figure |I| there the four- 
dimensional scalar curvature diverges in this limit. 
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w 



i i^ g in which case we obtain 

\a\g-2 



1 - 



\a\g-2 



'\a\b+\a\ 



H / Hg 

b V |o|g-2 



1 + 



' Hg 


\/\a\b+\a\ 


\a\g-2 





r 
r 



(6.6) 



N 
iVo 



1 + 



\ a \g 

\a\g-2 



Hg 
W\g-2 



2VFT 



(6.7) 



As seen in figure |] this solution contains three branches corresponding to physically different 
solutions. 



The right branch has g > 4r and r varies from to oo. In particular, g(r = 0) = oo and 
g(r = oo) = that corresponds to / = 1. It corresponds to — < W < ~ 7ypj^ - As # > |af we 
will soon see that this branch corresponds to negative masses. 



The center branch has < g < and r varies from j^ro to oo. In particular, g = ^ 
that corresponds to f (r — j^| r o) = and g(r = oo) = ^ that corresponds to / = 1. It 



2_ 

\a\ 



corresponds to — oo < W < — A. As g < j^- this branch corresponds to positive masses and 
is therefore interesting. Moreover, as can be seen from figure [I] for this solution the scalar 
curvature is everywhere regular. 



The left branch has g < and r varies from to ro- In particular, g (0) = oo and g(r = 0) = 0. 
It corresponds to ^= < W < 0. This branch has g < but no asymptotic region. It could 

y/\a\b 

however be patched together with the left branch on figure |3] to make a single geometry with 
< r < oo. 



The large distance expansion of the center branch is 



b 

\a\ 



N = 1 



2b ( Vb 
\a\ + VW> V Vb 

2h (Vb- 



\ + W\a\b 



r 



+ 0{r~ 



\a\ \ l°' + vl»l» r 



a + \/\a 



i\b \Vb 



+ 0{r~ 



(6. 



(6.9) 



where we chose N = ( ^ 



1 / b 

2 V lal 



This is full agreement with fl4.1|) and 
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7. Study of the solutions for a = 



In this case w = or oi = 2£. This situation is at the boundary of the range allowed in |§. Indeed in 
this case the term quadratic in the spatial derivatives of the scalar mode has a vanishing coefficient. 
The three basic equations now read 



rN' r 2 N 12 r 2 N" 
1 - / - r/' - (4/ + r/')^- + - 2/l£- 

rN > r 2 N t2 

r AT' r 2 A^' 2 r 2 AT" 







The second equation is solved for 



rN' 
~N~ 



N 



-1 + 



AT2 



N 



e 2 = l 



(7.1) 

(7.2) 
(7.3) 

(7.4) 



The other two equations are identically satisfied if ( |7.4[ ) is satisfied. Therefore the a = case is a 
degenerate case where f(r) is arbitrary while A" is given by (|7.4| ) . 

This reflects the degeneracy in the quadratic term for the scalar mode. 



8. Non-zero cosmological constant 

We will now proceed to understand solutions when the bare cosmological constant is non-zero. The 
relevant action now is 



S = J dtd 3 x^gN[a(K i:j K ij -\K 2 ) + f3CijC ij +£R + 01(040*) + a] 
The independent equations are 

+(w + 2£ > )(2w£ + 2(w + Oar 2 - wr£f') 2 = 



and 



rN' r 2 N' 2 
4£(/ - 1) - 2ar 2 + 8£/— + 2(w + 2^)/-^- = 



Once ( |3.2| ) is solved for / then N can be determined from ( |8.3| ) as 



rN' b 
~W ~ ~2 



where we defined 



2w 
(w + 20 



1 -4/1 



4£ 
(w + 20 



1 c 



/ 2b 



4a 
(w + 20 



a + b = 2 



(8.1) 



•2) 



(8.3) 



•4) 



(8.5) 
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As before, it suffices to consider one of the signs in front of the square root. 



To solve (|8.2|) we define 
so that equation ( $.2\) becomes 



f = ! + - + — r 2 + g(r) 
a ab 



r 2 g' 2 + 2a(rg' + g)g+(2b + cr 2 )g = 



(8.6) 



(8.7) 



One obvious generic solution is g = 0. Note that a can be completely rescaled out of this equation 
by redefining the radial variable as r 2 — > — . The only thing that remains in such a case is the sign 
of a. We will assume it to be positive (AdS), although similar remarks apply in the dS case. 
As before this equation can be simplified by the substitution 



to become 



(2b + cr 2 )W 2 rN' 
1 + abW 2 ' ~N~ 



2rW + (1 + abW 2 ) ^1 



'2a 2 b 3 W 2 + cr 2 (cr 2 + b 2 ) 



b 2 (cr 2 + 2b) 



(a + 2)cr 2 + 2afe T 



W = 



5.8) 



.9) 



cr 2 + 2b 

There is also another branch that is obtained by W — > —W that does not affect g and can be therefore 
ignored. Unfortunately, we have not managed to solve this equation exactly. It is however amenable 
to numerical integration straightforwardly. 

We will study however here its large distance asymptotics that will provide some interesting clues 
on the physics in the presence of a vacuum energy. This study is not trivial as the equation in question 
is non-linear and typically several long-distance asymptotics exist. 

8.1 The large distance asymptotics 

To find the large distance asymptotics, we will concentrate at distances r 2 3> | ^ | . In this region we 
can neglect the 2b term in equation ( |B.7| ) to obtain 

r 2 g' 2 + 2a{rg' + g)g + cr 2 g = (8.10) 

This equation can be solved exactly by the transformation 

cr 2 Y 2 



9 



1 + abY 2 
to obtain the equivalent equation 

2rW + (1 + abW 2 ) (1 - (a + 2)W) = 



(8.11) 



(8.12) 



which is the large r limit of ( |8.9| ). This can be solved exactly, and we will focus in the physical branch, 
a < 0. The general solution is 
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a + a_ 

0.38 r , i 




Figure 5: Left: a + as a function of \a\. Right: a_ as a function of \a\. 
Integrating for a < we obtain 



r 
r 



1 + ^\bY 



|a|-2)y + l|^T (8.13) 



with 



y/\a 


b±(\a 


-2) 


6 


a 


-4 



a± = "^T^M 1 ( 8 - 14 ) 



The values of a± for negative values of a are plotted in figure (|5|). We note that a + is always 
positive while a_ is negative for < \a\ < ~. 



From (p. 13) we observe that r — > oo on two distinct occasions. 



1. Y" — > i — and a_ > 0. This happens when lol > f . Expanding the solution we obtain 

'\a\b ' " 



Y=-^= + 0(r «-) , s~r 2+ - (8.15) 



2. V — > 2zi^r an d l a l > 2 or |a| < |. Here the asymptotics are 

, (8.16) 

Therefore 

• for < \a\ < |, the only possible large- distance asymptotics is g ~ r 2 with subleading terms 

3|a|-2 

that are powers of r h- 2 . 

• for | < |a| < 2 the only large- distance asymptotics possible is g ~ r 



finally for |a| > 2 both of the previous two asymptotics are possible. 
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The case a — — | must be treated separately. For this value the solution ( |3.13|) becomes 



= g 4(4^-3) 

»"0 



4F + 3 



4F-3 



(8.17) 



To obtain r — >■ oo we must take K — > | and we obtain asymptotically 



Y->-- 



9 



3c 2, 

H , £ ->■ lo S r + ■ 



(8.18) 



4 16 log r 

With the knowledge of the complete lists of asymptotics we can now proceed to establish the 
large distance expansions. 

8.2 The large distance expansions 

We start from solutions that asymptote to r 2 at large-distance. The general large distance asymptotics 
is as follows. 4 Solutions to ( |8.7| ) (with c=l) have the following large-distance structure 



119) 



n=0 



where go is given by 
9o(r) 

and in general 



+ 1_ _ 4(2 + 3a) 1 _ 8(a-2)(2 + 3a) 2 1 + 



6a + 4 r 2 (a + 6) r 4 (6 + a)(10 + 3a) r 
9n(r) = C n 



n 

1 -)- ^ ^ ' 



,2m 



m,=l 



In particular, the first few terms are 

2(2 + 3a) 2 1 2(-2 + a)(2 + 3a) 2 (4 + 28a + 25a 2 + 4a 3 ) 1 



9i(r) = C 



1 - 



(2 + a) 2 r 2 

(2 + 3a) 2 



92{r) 



(2 + a) 



C 2 



(2 + a) 4 (6 + a) r 4 

:| 2a(a-2)(2 + 3a)(10 + 7a) 1 
+ (2 + a) 2 (6 + 5a) ^ + 



+ 0{r- 



(8.20) 
(8.21) 

(8.22) 
(8.23) 



4a(a - 2) 2 (2 + 3a)(768 + 4160a + 7288a 2 + 5500a 3 + 1786a 4 + 185a 5 ) 1 mt 

H — \T7^ — . Z ^777. . - 7 + 0{r 



(2 + a) 4 (6 + a)(6 + 5a)(10 + 7a)) r 4 

C is an integration constant and C n ~ C n in ( |8.21| ). This expansion is well defined for generic 
irrational a. For rational values of a resonance effects can happen and logarithms may appear in this 
expansion. 

From (|3.4|) we may compute the large distance expansion of the lapse. We find 

- 2 a(a-2), (-2 + a) 2 a(-32 -44a + 4a 2 + a 3 ) 



log N = h (a - 2)r + 



logr 



16(2 + a) 2 r 2 



(8.24) 



4 We rescale r so that we set c = 1. 
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(a-2)V(3a + 2) 2 C _», a-a 

-7- ^a + a _)_ 



(a + 2)3(a + 4) 

For the blackness function /, we observe that this class of large distance asymptotics have the 
standard r 2 behavior characteristic of gravity with a cosmological constant, albeit the coefficient of 
the r 2 term is affected by the new term in the action. Indeed in terms of the original couplings £, a\ 
and a, the blackness function / from (|8.6|) has an expansion 



24 - a\ ? (34 - 2oi) (24 - ai) 

that reduces to the standard asymptotics when a x = (corresponding to a — > —00). 

What is remarkable in (|8.19|) is that no 1/r tail exists generically. Indeed the leading term 



proportional to the integration constant behaves as r 5 ^, and only becomes 1/r in the GR limit 
a — > —00. 

The lapse on the other hand has completely different large distance asymptotics that are expo- 
nential to leading order, and where the "mass" contribution (proportional to C in ( |B.24 ) is strongly 
subleading. 

It is therefore expected that the interplay between a small cosmological constant and the large 
distance tails of spherically gravitating objects will be non-trivial, and unfortunately is not easily 
accessible through a large distance expansion alone. On the other hand, it is expected that there may 
be non-trivial constraints arising from this interplay if the cosmological in the universe is assumed to 
be non-zero. 

We now move to the other alternative large distance asymptotics described in the previous sub- 
section. In this case the solution can be written as 

S(r) = i> 2+ " 9n(r) , 9n(r) = C n (l + J2^A (8-26) 

n=l \ m=l / 

C\ is the integration constant while the other coefficients are non-linearly related to it. The lapse 
has the same leading behavior as in ( 8.24J) . Note that in the phenomenologically interesting region 



a — > —00, l/a_ is large and positive. Therefore this branch of solutions have a large power dependence 
which make it phenomenologically uninteresting. 



Notes Added 



• While this work was written, || has argued that the modified Hofava-Lifshitz gravity proposed 
in H has also potential strong coupling regions, based on the two derivative approximation. 
These regions are not visible in the spherically symmetric solutions studied here. Subsequently 
in |0] it was argued that the breakdown of the two derivative theory is prevented by the higher 
derivative terms rendering a theory without strong coupling problems. 
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It was also brought to our attention after the completion of this paper that the equations solved 
for zero cosmological constant are identical with those governing a special class of solutions 
of a Einstein- Aether theory [jnj, with c\ = C3 = 0, but with the aether vector restricted to 
be hypersurface orthogonal. This is the case with spherically symmetric solutions. Indeed the 
solutions found here with a < match those found in ||, and their linearized versions studied 
earlier in ||. It was also observed that the physical solution corresponding to the center branch 
of figure [| can be extended inside the horizon by the right branch of the solution in figure [31 In 
that case the horizon at / = is a minimal sphere. The relation with Einstein- Aether theory 
was also discussed in [[F|. 

After the appearance of this paper, it was argued in []12 |, echoing an earlier claim in [|TTJ that 



the second class constraints are very strong and in generic solutions they force the lapse to 
vanish. This does not happen in solutions with symmetry like those that are treated here, but 
this is certainly an important issue that deserves further attention. 



After the completion of this work we became aware of reference |13| where spherically symmetric 



solutions in standard Hofava-Lifshitz gravity with non-zero shift were presented. 
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APPENDIX 



A. Study of the solutions for a > 



We present in this appendix the solutions to the equations without a cosmo logical constant for a > 0. 

In this branch, as a > it implies that b < 2, and therefore corresponds to the linearly unstable 
class, 0. 

The solution is 



y/l + abW 2 exp 


yb/a 


aictan(y/abW) 




1-aW 





r 
r 



while for N we obtain 



N = N exp 



—\ -arctan(v abW) 



The solution contains two branches that are visible in figure [5| 



N 
N 



15 




10 




5 





-3 -2 -1 



(a.i; 



(A.2) 




Figure 6: On the left r/ro versus W for a = b = 1. On the right N/Nq versus W for a = b = 1. 

Substituting for g we obtain 
For IF = J^J For this to be real > ag > -2. 



2+ag 



exp 



arctan 



2+ag 



-ag 



r 



N = N exp 



b V 2+ag 

The functions g(r) and N(r) are shown in figure |7|. 



-arctan 



-ag 



2 + ag 



(A.3) 



The left branch — - < q < — - has an asymptotic infinity at q = — - and since it remains 

n. n. j sr j ij ft 



smaller than this everywhere it has a positive mass parameter. It corresponds to 
W < oo. 



< 
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Figure 7 : On the leftr/ro versus g, with W > for a = b = 1. On the right, N as a function of g normalized 
to one at r — > oo. 

— The right branch — - < g < has an asymptotic infinity at g = — - and since it remains 
larger than this everywhere it has a negative mass parameter. It corresponds to < W < 



For the large distance limit g — > — - we obtain 



b , 

h o\j -exp 

a V a 



arctan 



r 



+ 0{r~ 



(A.4) 



N 2 



1 + b\ — exp 

a 



— arctan 
a 



^ + 0(r 
r 



(A.5) 



where we chose N = exp 



b arctan ' ' " 



For W 



2+ag ' 



For this to be real > ag > —2. 



2+ag 



exp 



- arctan 

a 



-g-g 

2+ag 



1 + 



2+ag 



r 



N = N exp 



oil -ag 

arctan \ / 

a V V 2 + ag 



(A.6) 



Here as shown in figure |8| there is no asymptotic infinity. The radius is bounded < r < r$. It 
is quite plausible that this solution combines with the right branch of figure [7|, as suggested by 
figure H in order to make a single solution valid from < r < oo. This solution has a negative 
mass parameter. 
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g 



Figure 8: On the left r/r® versus g, with W < for a = b = 1. There is no asymptotic infinity and r varies 
from to Tq. On the right, N/Nq as a function of g. 
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